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The electron, hole, and exciton ^-factors and diamagnetic coefficients have been calculated us- 
ing envelope- function theory for cylindrical InAs/InP quantum dots in the presence of a magnetic 
field parallel to the dot symmetry axis. A clear connection is established between the electron g- 
factor and the amplitude of the those valence-state envelope functions which possess non-zero orbital 
momentum associated with the envelope function. The dependence of the exciton diamagnetic coef- 
ficients on the quantum dot height is found to correlate with the energy dependence of the effective 
mass. Calculated exciton ^-factor and diamagnetic coefficients, constructed from the values associ- 
ated with the electron and hole constituents of the exciton, match experimental data well, however 
including the Coulomb interaction between the electron and hole states improves the agreement. 
Remote-band contributions to the valence-band electronic structure, included perturbatively, reduce 
the agreement between theory and experiment. 



I. INTRODUCTION 



The magnetic moment and Zeeman energy splittings of 
an electronic spin are controlled by the g-tensor, whose 
value differs from the bare electron's ^-factor of ~ 2 due 
to the spin-orbit interaction in a solidi. Thus an under- 
standing of the g-tensor not only provides insight into 
the spin-orbit interaction in a solid, but it also can yield 
insight into the effective orbital motion (and resulting an- 
gular momentum) associated with a quantum state 1 . The 
spin-orbit interaction affects the magnetic moment of an 
electronic state by permitting a contribution from the 
orbital angular momentum in addition to the spin mo- 
mentum (leading to an effective Lande ^-factor in atomic 
systems). In solids, in addition to the atomistic (Bloch) 
orbital angular momentum, the envelope wave function 
also contributes to the magnetic moment via the envelope 
orbital angular momentum. 

The ^-factors of bulk semiconductors, which often dif- 
fer substantially from the electron's bare ^-factor, can- 
not be explained solely by the large orbital angular mo- 
menta of Bloch states that contribute to the orbital mo- 
tion of the state within a crystal lattice unit cell. In- 
stead the spatial extension of the wave function across 
many unit cells leads to the large ^-factors through in- 
teratomic circulating currents originating from the spin- 
orbit potential 2 -. As a wave packet in a semiconductor 
crystal can have a substantial extent due to the small 
effective mass m*, the envelope wave function can sus- 
tain large extended circulating currents leading to a large 
envelope orbital angular momentum. In the absence of 
spin-orbit coupling these currents vanish and do not con- 
tribute to the modification of the ^-factor at all. How- 
ever, in most semiconductors the spin-orbit coupling is 
substantial and will lift the spin degeneracy, leading to 
a finite spin-dependent envelope orbital angular momen- 
tum. Therefore the ^-factor is modified from the Lande- 



factor by means of the additional envelope orbital angu- 
lar momentum generated by circulating currents in the 
envelope wave function. 

These circulating currents are reminiscent of the ori- 
gin of diamagnetism: the response of electrons to an 
applied magnetic field is to create circulating currents, 
leading to a magnetization opposed to the applied field. 
The underlying origin of the diamagnetic circulating cur- 
rents is different than those associated with the g-tensor, 
which come from the spatially-periodic spin-orbit poten- 
tial. However, for both phenomena the generation of 
orbital angular momentum by extended circulating cur- 
rents is essential to understand the effect and explain the 
experimental results. It was Ehrenfest who initially iden- 
tified the role of extended circulating currents in produc- 
ing large differences in the diamagnetic susceptibilities of 
various materials^. 

In quantum dots, confinement alters the envelope or- 
bital momentum of the envelope wave function by pre- 
venting it from extending over as many lattice unit cells 
as it would in the bulk. In the limit of infinite confine- 
ment potentials and vanishing dot size, the envelope or- 
bital momentum is quenched and the ^-factor is solely 
determined by the atomistic (Bloch) Lande-factor. This 
idea has been successfully applied to paramagnetic im- 
purities in salts^. The effect of envelope orbital momen- 
tum quenching was demonstrated in Ref. HI for a spher- 
ical, unstrained InAs nanocrystal, by showing that the 
electron ^-factor approached the bare electron value of 
~ 2 much more quickly than would be expected from 
a parametrized bulk model based on Ref. 6, A simi- 
lar quenching effect on the orbital angular momentum 
should reduce the diamagnetic coefficient in small quan- 
tum dots with large barriers. For both quantities it is 
crucial to explore how confinement affects the envelope 
orbital angular momentum. An understanding of the 
nature of the ^-tensor's dependence on dot composition 
and structure can thus clarify the role of orbital angu- 
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lar momentum quenching in determining this property 
in a quantum dot 5 . Diamagnetic coefficients provide a 
picture of the orbital motion associated with a quantum 
state through a different path than the g- tensor; a the- 
ory which correctly describes both has passed a stringent 
test of its picture of the orbital motion of the electronic 
state. 

In addition to these fundamental concerns, a single 
spin in a quantum dot provides a model system for the 
observation and manipulation of an individual quantum 
system^, and has been proposed as a qubit for use 
in quantum computation^, ^-tensors (but not diamag- 
netic coefficients) figure prominently in spin-based pro- 
posals for quantum information processing as they pro- 
vide an effective spin-manipulation method, and hence 
qubit-operation mechanism. For example, proposals for 
single-spin control via an electric field include electric 
field modification of the g-tensor to bring the spin into 
and out of resonance with a global RF field^"— , or to 
tilt the spin's precessional axis^ 2 -"— . They also provide 
energy-level structures favorable for optical pumping of 
spin polarization of an electron^—, or a hole— ^ in a 
quantum dot. Due to energy-conservation constraints 
the g-tensor also can control the spin lifetime (Ti) for 
a dot* 2 ^. 

Compositionally-defined quantum dots 24 permit much 
larger barrier heights than gate-defined dots, as well 
as efficient optical access with simultaneous control of 
charging 25 . The larger confinement energies in such dots 
for both electrons and holes from these barriers may in 
principle permit single-spin behavior at higher operat- 
ing temperature than in lithographic dots^, however the 
intrinsic strain and dot shape asymmetry provide a chal- 
lenge to theory. Shown in Fig. [I] is a schematic dot, 
with height h and radius r, which can be formed through 
metal-organic vapor-phase epitaxy (MOVPE) of InAs on 
an InP substrate. Details of the growth of these roughly 
cylindrical InAs dots can be found in Ref. [26|. Recent ex- 
perimental work^i has investigated the magneto-optical 
properties of such dots by probing the component of the 
exciton g-tensor along the [001] -direction (this compo- 
nent is from now on referred as the ^-factor). A strong 
dependence of the exciton ^-factor (g ex ) on the dot height 
was shown, whereas the same work found the exciton dia- 
magnetic coefficient (a ex ) was mostly independent of the 
height. The quantum dot radii had little influence on 
either g ex or a ex . Theoretical calculations are needed to 
provide insight into the details of the observed depen- 
dences of g ex and a ex and to provide clarity about the 
ability to tune the g ear factor through control of the size 
and shape of the quantum dot. 

Here we describe an extensive series of multi-band 
envelope-function calculations of electron, hole, and ex- 
citon ^-factors and diamagnetic coefficients based on an 
eight-band strain-dependent k • p model of the bulk semi- 
conductor constituents for a series of quantum dots sim- 
ilar to the ones reported in Ref. l27|. The calculational 
method is detailed in Sec. Hfl The calculated ^-factors 




FIG. 1. (Color online) The quantum dots are made out of pure 
InAs (red) embedded in InP (grey), and have a cylindrical 
shape with a radius r and height h. 



and diamagnetic coefficients of the electron and hole are 
discussed separately in sections IIIII and IIV1 The exciton 
^ex-factor and a exi which are constructed from the values 
for the electron and hole, are compared to experimental 
data in section [VI The Appendices provide a detailed 
comparative analysis with our results of other schemes 
to calculate the ^-factor and diamagnetic coefficient. 

II. METHODS 

The ^-factors and diamagnetic coefficients of quantum 
dots are calculated using eight-band strain-dependent 
k • p-theory on a real space grid 2 ^. The strain is cal- 
culated using linear elasticity continuum theory. All cal- 
culations are performed at T = K and material pa- 
rameters are taken from Ref. 29. The magnetic field is 
coupled both to the envelope wave function and the spin 
of the Bloch functions 5 , where the latter is taken into 
account by the Zeeman Hamiltonian: 

1 /2<r \ 

^Zeeman = ~^ B B • §J , (1) 

V w 

where \±b is the Bohr magneton and a and J are the spin 
matrices of spins | and | respectively. The ^-factors for 
the conduction, valence and spin-orbit bands are respec- 
tively 2, | and | (i.e. the Lande- factors). Note that we 
have explicitly left out the remote band contributions. 
This approximation is validated in appendix lAl 

The calculated quantum dots have a cylindrical shape 
and are assumed to be pure InAs (see Fig. [I]), as 
was inferred from cross-sectional scanning tunneling mi- 
croscopy in Ref. H3. Calculations for various heights h 
and radii r for the quantum dots have been performed 
to study the size dependence of the ^-factor and diamag- 
netic coefficient. The energy levels and wave functions of 
these dots are calculated with magnetic fields of — 10 T 
applied in the growth [001] -direction. The sizes of the 
quantum dots are small enough compared to the magni- 
tude of the applied magnetic field to preclude any transi- 
tion to complete magnetic confinement, which has been 
studied elsewhere^. The spin-dependent energies of both 
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FIG. 2. (Color online) Height dependence of the electron g c -factor of (a) unstrained and (c) strained InAs/InP quantum dots, 
respectively, (b) the correlation between the g c -factor and the sum of all components having L e z nv = dzl for unstrained QDs. 
(d) the height dependence of the sum of all components having envelope angular momenta L e z nv = ±1 (dots) and L e z nv = 
(squares). The different colors indicate different radii of QDs. Horizontal lines are the g c -factors of quantum wires having 
different radii. 



the electron (c) and hole states (v) can be parametrized 
as 

E c (B) = E c , ± X -\vb\9cB + a c B\ (2) 

E v (B) = E Vj0 =F \\ii B \ {g° v + g 2 v B 2 ) B + a v B\ (3) 

where E^ v ^q is the ground state energy, g( c ,v) the g- 
factor, fiB the Bohr magneton, and ot( CjV ) the diamag- 
netic coefficient. Note that the hole (but not the elec- 
tron) Zeeman energy is observed to possess significant 
nonlinearity; the hole ^-factor has a zero-order compo- 
nent g® and a second-order component g%. The upper 
(lower) signs in front of the Zeeman energy correspond 
to spin states parallel (anti-parallel) to the applied mag- 
netic field. Using these definitions for the energy levels, 



we can define the magneto-optical properties of the exci- 
ton as: 

9ex = -Qc- 9v, (4) 

&ex = \0L C \ + \ot v \, (5) 

E em = \E Cj o\ + |^,o | , (6) 

where g ex is the exciton ^-factor, a ex the exciton dia- 
magnetic coefficient, and E em the emission energy of the 
QD. These definitions are consistent with the experimen- 
tal definitions of Ref. H3- 

III. ^-FACTORS 

The calculated electron g c -factors and hole ^-factors 
are shown in Figs. [2] and [3j The different colors indicate 
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different radii r of the QDs. The continuous horizontal 
lines are the ^-factors calculated for quantum wires (i.e. 
quantum dots having infinite height). 

A. Electron state 

For the electron state, we have computed two cases: 
without strain, Fig. E](a), and with strain, Fig. [2](c). In 
both cases the g c -factor has a similar trend both in terms 
of height and radius; as confinement increases by decreas- 
ing either height or radius, the ^-factor tends towards the 
free electron ^-factor value of +2 due to orbital angular 
momentum quenching. As explained in section [TJ the 
envelope orbital momentum causes the deviation of the 
electron ^-factor from the Lande-factor. It is therefore 
necessary to analyze the origin of the envelope orbital 
angular momentum in more detail. We first focus on the 
results for unstrained quantum dots. 

Within the framework of eight-band k • p-theory, any 
quantum state consists of four different band compo- 
nents: conduction band (CB), heavy- hole (HH), light- 
hole (LH), and split-off (SO). Each of these components 
has a Bloch and envelope wave function with correspond- 
ing (total) orbital angular momenta J Bloch and L env . Due 
to strong spin-orbit coupling these momenta couple with 
each other giving rise to a combined total momentum F: 

\F, F z ) = | J Bloch , Jf och ) |L env , Lf v ) (7) 

The calculated electron ground state has an almost com- 
pletely F = | character, which can be formed from 
the four different band components by taking the proper 
combinations of J Bloch and L env . However, not all L env 
are allowed; the parity of the effective mass equations 
requires that only either even or odd L env can be mixed 
into the stated. All allowable combinations have been 
listed in table [TJ together with the probability density of 
the different components along the [001] -direction. 

Table U shows the CB component has a L env = and 
therefore does not contribute to the envelope orbital an- 
gular momentum and thus to the deviation of the g c - 
factor from +2. However, the valence band portions of 
the wave function have components with |L env ,L| nv ) = 
|1,±1), which do lead to a finite envelope orbital an- 
gular momentum directed along the [001] -direction. In 
Fig. EJb) we have plotted the g c -factor against the sum 
of all components of the wave function having L e z nv = ±1. 
For fixed radius, there is a nearly linear relation between 
the g c -factor and the amount of the wave function which 
has an envelope orbital momentum L e z nv = ±1. Less am- 
plitude in the band components having L e z nv = ±1 leads 
to a g c -factor tending to +2. In other words, the value of 
the electron g c -factor for a particular quantum dot hav- 
ing a certain height and radius is solely determined by 
how the confinement has affected the composition of the 
wave function. 

To investigate this further, we have plotted in Fig.[2jd) 
the sum of all components having either L e z nv = ±1 (dots) 



and L e z nv = (squares) as function of height. The ob- 
served dependencies can be straightforwardly understood 
by examining the underlying Hamiltonian 28 . The cou- 
pling between the CB-component and the valence band 
components having L e z nv = ±1 is proportional to k x and 
k y , whereas the coupling with the components having 
L e z nv = is proportional to k z . The allowed k- vectors for 
our cylindrical quantum dots having height h and radius 
r are roughly proportional to k z ~ l/h and k Xl k y ~ 1/r, 
and thus for larger height the fraction of L e z nv = de- 
creases. Although not directly linked to the height, the 
fraction of L e z nv = ±1 increases simultaneously. This is 
merely due to the relative increase of the 1/r coupling 
(L| nv = ±1) compared to the directly affected l/h cou- 
pling (L e z nv = 0). In a similar fashion the radial de- 
pendence can be understood. This analysis shows that 
quantum confinement is only affecting the composition of 
the wave function via the allowed k- vectors in the quan- 
tum dot. Thus the effect of quantum confinement can 
be understood as probing a bulk material at a specific k- 
vector determined by the size of the quantum dot. This is 
analogous to earlier work on quantum wells which found 
the non-parabolicity of the electronic effective mass by 
varying the quantum well thickness. 

As can be inferred by comparing the unstrained and 
strained ^-factors in Figs.^a) and (c), strain is changing 
the g c -factor only quantitatively and not qualitatively. 
The strain Hamiltonian has the same functional form as 
the kinetic part of the Hamiltonian 32 , meaning that its 
effect will only be able to modify the size of the couplings 
and energy gaps between the different bands. From this 
point of view strain can be regarded as a perturbation 
onto the system, not affecting the mechanisms determin- 
ing the g c -factor. In particular, compressive strain inside 
the QD will cause the effective band gap to become larger 
compared to unstrained QDs. This gives rise to an overall 
smaller fraction of the valence band mixed into the elec- 
tron state, and so for similar sizes, the strained g c -factors 
have values closer to +2 than the unstrained ones. 



B. Hole state 

The hole state presents additional complications versus 
the electron state. We show calculations without strain, 
Fig. [31(a), and with strain, Fig. EJc). As for the electron 
state, the envelope orbital angular momenta of the four 
different band components of the wave function are es- 
sential to understand the behavior of the ^-factor. As 
the hole ground state has a strong F z = d=| character, 
different envelope orbital momenta L env are mixed into 
the state. The proper combinations of J Bloch and L env 
are listed in table El together with probability density 
of the different components along the [001] -direction. As 
can be seen from the probability densities, higher orbital 
angular momenta are mixed into the hole state, since the 
energy gaps between the valence bands are rather small. 
There are a number of reasons why an analysis of the 
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TABLE I. Combinations of J Bloch and L env leading to \F, F z ) = |±, ±±), for each of the different components % of ^ (r), of a 
QD having height of 6 nm and radius of 11 nm. The color scale is different for each plot. 
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FIG. 3. (Color online) The height dependence of the zero-order hole gy-iactor of ( a ) unstrained and (c) strained InAs/InP 
quantum dots, (b) the dependence of the zero-order hole g£-factor on the LH-component for unstrained and strained InAs/InP 
quantum dots, (d) the height dependence of the second-order hole ^-factor of strained InAs/InP quantum dots. The different 
colors indicate different radii of QDs. The horizontal lines are the g£' 2 -factors of quantum wires having different radii. 



envelope orbital momenta is more complicated than for 
the electron state: 

No clear limit at infinite confinement. The HH- 
component has both L env = and L env = 2 components, 
each of which has L e z nv = and thus does not contribute 
to ^-factor. However, the L env = and L env = 2 parts 
of the wave function will have a different Lande-factor. 
Depending on the relative weight of these orbital angu- 
lar momenta, the ^-factor at infinite confinement (in the 
completely quenched limit) will vary. Thus we do not ob- 
serve the clear limit at small heights like was obtained for 
the electron state. 

Strain is no longer only perturbative. The electron g c - 
factor is determined by the admixture of valence bands 
into the electron state, a mechanism upon which strain 
only acts perturbatively. For holes the strain plays a 
more dominant role, and changes energy gaps and cou- 



plings to a much greater extent. Indeed by comparing the 
unstrained and strained ^-factors [Fig. [3^a) and (c)], it 
can be seen that even the qualitative trends are not the 
same. 

Smaller ground- excited state energy splitting. At 
heights larger than ~ 15 nm an excited state crosses the 
ground state, leading to mixing between the two. This 
means that also envelope orbital angular momenta con- 
tributions get mixed, complicating the analysis. 

Due to these complications a single driving mechanism 
for the ^-factor cannot be extracted. However, some 
general trends can still be observed. There is a very dif- 
ferent trend of the ^-factor as a function of height for 
different radii. At large heights there seems to be an 
overall converging trend. This converging trend is more 
dominant in unstrained quantum dots [Fig. 03(c)] . The di- 
vergence, more pronounced for large radii quantum dots, 
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is enhanced by strain [Fig. 09(b)]. 

Besides the usual zero-order ^-factor, we have ob- 
served that there is a considerable second-order ^-factor 
affecting the Zeeman energy of the hole states [Fig. 09(d)]. 
Note that this quantity has per definition a dimension 
(see Eq.[3]) and is thus dependent on the magnitude of the 
applied magnetic field. For example a QD with a height 
of 10 nm and radius of 13 nm has a g% ~ 6 • 10 _3 /T 2 lead- 
ing to a deviation of the zero-order ^-factor from ~ 4 to 
~ 4.6 at a magnetic field of 10 T. Such deviations are cer- 
tainly large enough to be noticed in experiments. To the 
best of our knowledge, however, such non-linearities have 
only been seen in very tall dots^. In quantum wells, non- 
linearities of the hole Zeeman energy have been observed 
experimentally and were reproduced theoretically-^—. 
In these studies the non-linearity was attributed to the 
admixture between HH- and LH-subbands. In that case 
mass reversal is induced by quantum confinement, lead- 
ing to an anti-crossing of the HH- and LH-subbands as 
function of magnetic field. This anti-crossing only oc- 
curs for one of the two spin split subbands, leading to a 
non-linear Zeeman energy. However, in our calculations 
both spin split ground states exhibit a non-linear depen- 
dence as function of magnetic field. Therefore the mech- 
anism as explained for quantum wells cannot explain the 
origin of the ^-factor in quantum dots. We attribute 
the non-linearity to mixing between the HH-like ground 
and LH-like excited states, and have verified that the 
magnitude of g% decreases with increasing strain (larger 
energy splitting) and decreasing 72 and 73 parameters 
(smaller coupling), both resulting in a smaller degree of 
mixing. Furthermore, we observe that the non-linearity 
has a strong size dependence: increasing either height or 
radius decreases the energy splitting between the ground 
and excited states, leading to a larger ^-factor. It re- 
mains unclear why this effect is strongly present in the 
calculation and has not been observed experimentally. 



IV. DIAMAGNETIC COEFFICIENTS 

We will now discuss the diamagnetic coefficient for the 
electron and hole state (see Fig. [4|). Within the frame- 
work of the k • p-approximation, in first order perturba- 
tion the diamagnetic coefficient of a carrier in a semicon- 
ductor is proportional to2i: 

a ~ ^ (8) 

where (r 2 ) is the average lateral extension of the wave 
function perpendicular to the applied magnetic field (and 
can therefore be directly associated with the radius), and 
m* the effective mass of the carrier. As a result there is a 
strong radial dependence of the diamagnetic coefficient, 
both for the electron and hole state: a large radius quan- 
tum dot has a large diamagnetic coefficient. The strong 
height dependence is unexpected; for the electron state 
the diamagnetic coefficient becomes smaller at smaller 



heights, whereas for the hole state it shifts both to larger 
and smaller values depending on the height. We will start 
by discussing the electron state. 

The height dependence should arise either from details 
in the average lateral extension (r 2 ) or the effective mass 
m*, as can be inferred from Eq. (|5J). It has been ver- 
ified that (r 2 ) is only weakly dependent on the height 
(< 10% for the most affected 15 nm radius quantum 
dots), and hence cannot be responsible for the strong 
deviations in a c of almost a factor of 2. We therefore 
attribute the strong height dependence to the energy- 
dependence of the effective mass. As mentioned in previ- 
ous sections, quantum confinement is effectively probing 
the bulk dispersion at specific fc-vectors. At the small 
quantum dot heights considered here, the dispersion is 
sampled at large k- vectors, where the conduction band is 
no longer parabolic and corrections to the effective mass 
have to be made, i.e. m* = m*(h). For InAs these cor- 
rections lead to a more heavier mass at large k- vectors, 
causing the diamagnetic coefficient to become smaller at 
smaller heights, as can be observed in Fig. [U We want 
to note that we have verified that strain is not an impor- 
tant mechanism in determining a c , as was the case for 
the electron g c -factor. The sole effect of strain is to mod- 
ify the magnitude of the diamagnetic coefficient (mainly 
via m*), the qualitative trends stay the same. 

The hole state diamagnetic coefficient has an even 
more peculiar height dependence: at intermediate heights 
it decreases, whereas at large heights it increases. The 
nature of the hole state induces a number of ways by 
which the diamagnetic coefficient can be modified. First 
of all the lateral extension can be influenced by size de- 
pendent HH-LH mixing. We have verified that (r 2 ) varies 
by < 15% for the 15 nm radius quantum dots (and varies 
less for smaller radii). Contrary to the electron state, this 
variation is of the same order as the deviations in a v . 
However, the effective mass will also vary with height. 
It is a priori not possible to disentangle these effects 
and study their respective influence. Note that we also 
calculated the diamagnetic coefficients including the k- 
parameter. Similar to the ^-factors, to first order the k- 
parameter does not modify the diamagnetic coefficients 
(see appendix [A] for a discussion about the remote band 
parameters). 



V. COMPARISON WITH EXPERIMENT 

Using the separate electron and hole ^-factors and dia- 
magnetic coefficients, one can compute the exciton g ex - 
factor and exciton diamagnetic coefficient a ex , as defined 
by Eqs. (|4])-(j6]). Since in the experiment the Zeeman en- 
ergy was found to be linear with magnetic field, we have 
only taken the g^-factor into account, which is exactly 
valid for small magnetic fields. Since experimentally only 
the emission energy E em is accessible, we have plotted the 
^e^-factor and a ex as function of E em in Fig. [5l Again 
the different colors represent different radii QDs. Note 
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FIG. 4. (Color online) The height dependence of the (a) electron diamagnetic coefficient a c and (b) hole diamagnetic coefficient 
a v of strained InAs/InP quantum dots, respectively. The different colors indicate different radii of QDs. The horizontal lines 
are the a's of quantum wires of different radii. 



that the different dots along a line of fixed radius repre- 
sent QDs with roughly a ML height difference. The black 
squares are the experimental data points from Ref. [27]. 
We will first discuss the g ear factor. 

It was experimentally established 27 that the general 
trend of the g ex -fdiCtoT with emission energy originates 
from the dependence on the height of the QDs. This 
is confirmed by the theoretical calculation, in which ir- 
respective of the exact radius of the QDs the height 
dependence matches with the experimentally observed 
trend. Furthermore, from cross-sectional scanning tun- 
neling microscopy analysis it was inferred that the max- 
imum radius of the QDs was ~ 15 nm, meaning that the 
13 nm radius theoretical calculations that overlap the 
experimental data are a good match. Experimentally 
the spread around the trend was related to variations in 
the lateral size: QDs below the trend had on average 
a smaller a ex than QDs above the trend. The oppo- 
site trend is observed in the calculations. However, the 
assumptions made about the shape and composition of 
the quantum dot may be too simplistic to reproduce this 
detailed trend; for example, a size-dependent composi- 
tion profile might alleviate the discrepancy. Note that 
it is mainly the hole ^-factor which determines the g ex - 
factor; the electron g c -factor is nearly independent of the 
radius when plotted as a function of E em . We want to 
stress that, since the g e;E -factor is a combination of the 
# c -factor and ^-factor, it is difficult to extract relevant 
information from it to identify the physics determining 
the ^-factors of the individual carriers. It is far more 
useful to have experimental access to the g c -factor and 
(^-factor separately, preferably in a system where the size 
of the quantum dot can be highly controlled. 

The diamagnetic coefficient a ex is plotted against E em 



in Fig. [He). Experimentally it was established that the 
emission energy is dominated by the height dependence. 
Assuming that a ex is then governed by the radius (see 
Eq. [8|), there should only be a weak dependence between 
a ex and E em . This has been observed experimentally, 
and confirmed the interpretation of the height-dominated 
E em . The theoretical calculations support this analysis; 
there is only a rather weak dependence between a ex and 
E ern . The calculations, however, show that the electron 
and hole diamagnetic coefficients separately are strongly 
dependent on the height (Fig. [4]). It is the coincidental 
height dependence of the emission energy and the (par- 
tially) opposing height dependence of the electron and 
hole diamagnetic coefficients, which conceal the true size 
dependency. Experimental access to the separate elec- 
tron and hole energy levels would be the only way to 
verify this. The match in terms of absolute size of the 
quantum dot is less satisfactory than for the ^-factor; the 
calculations with a radius of ~ 8 nm fit the diamagnetic 
coefficient measurements best. 

To improve on the match for a exi we included the 
Coulomb interaction between the electron and hole states 
via a self-consistent iterative Hartree approximation^. 
After two iterations more than 99% of the total modifi- 
cation to the emission energy is taken into account. In 
graphs (b) and (d) of Fig. [5] the influence of the Coulomb 
interaction to respectively g ex and a ex are displayed. Ex- 
cept for corrections to the emission energy (< 7%), there 
are only small modifications to g ex . This was not too 
surprising, as the Coulomb interaction might modify the 
spatial extension of the wave function, but not as much 
the composition of the states. Supporting that picture, 
the diamagnetic coefficient is much more strongly modi- 
fied: it becomes smaller, such that there is a better agree- 
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FIG. 5. (Color online) The emission energy dependence of the (a) exciton g ex -factor and (c) diamagnetic coefficient a ex of 
InAs/InP quantum dots, respectively. The emission energy dependence of (b) g ex and (d) a ex with Coulomb interaction taken 
into account, respectively. The different colors indicate different radii of QDs. 



merit between the experimentally found radius (~ 13 nm) 
and the theoretical match (~ 11 nm). The Coulomb 
interaction pulls the electron and hole wave functions 
together, such that the exciton state has a smaller lat- 
eral extension (and therefore a ex ) compared to the non- 
interacting states for a similar radius quantum dot. 

VI. CONCLUSIONS 

We have calculated the ^-factors and diamagnetic coef- 
ficients of electron and hole states in InAs/InP quantum 
dots. The electron g c -factor is determined by how quan- 
tum confinement affects the amount of envelope orbital 
angular momentum mixed in through the admixture of 
valence-band contributions. In the limit of infinite bar- 
riers and a dot of vanishing size the electron g c -factor 
will approach the Lande- factor of +2. Strain acts merely 



a perturbative mechanism on the g c -factor. A similar 
analysis is less straightforward for the hole ^-factor. We 
attribute the observed non-linear hole Zeeman energy to 
ground-excited state mixing. The electron and hole dia- 
magnetic coefficients are determined both by the lateral 
extension of the wave function and the effective mass. 
The height dependence of the electron diamagnetic coef- 
ficient is due to a strongly height- (or energy-) dependent 
effective mass. Using the single-particle states, the calcu- 
lated exciton ^-factor and diamagnetic coefficient agree 
well with the experimental data. Including the Coulomb 
interaction between the electron and hole states improves 
the agreement, especially by decreasing the diamagnetic 
coefficient due to reductions in the lateral extent of the 
wave function. 
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FIG. 6. (Color online) Height dependence of the zero-order hole g v -factor of InAs/InP quantum dots (a) without and (b) with 
the ^-parameter included, respectively. Different colors indicate different radii of QDs. The horizontal lines are the ^-factors 
of quantum wires of different radii. 



Appendix A: Remote band contributions 

A k • p-calculation takes only a finite number of bands 
into account, and the influence of remote bands can be 
incorporated via phenomenological parameters. These 
remote band parameters then do not change as the con- 
finement experienced by an electron or hole in the dot 
changes. As explained in section H orbital angular mo- 
mentum quenching strongly influences the ^-factor and 
diamagnetic coefficient. The influence of remote bands, 
considered perturbatively in a k • p-calculation, is not 
quenched with increasing confinement. The quantum 
dots considered in this work are in the strongly confined 
limit, so it is reasonable to assume that the contribu- 
tion of the remote bands has been completely quenched. 
In the calculations presented above we have left out the 
remote band parameters for the ^-factors and used the 
Lande-factors in the Hamiltonian. This is certainly a 
very good approximation for the conduction band g- 
factor, since in In As remote bands contribute only about 
4% to the bulk electron g c -factor— . For holes the situa- 
tion is more complicated, since the remote band contribu- 
tion is non-negligible and one should consider including 
the phenomenological ^-parameter— in a bulk calcula- 
tion. 

To view the effect of ft, the calculations are performed 
both with and without it [Figs. [6](a, b)]. The dependence 
of the ^-factor with height is not very different for both 
cases, however including the ^-parameter means that a 
large offset is introduced (the unquenched effect of these 
remote bands). This offset is very large compared to 
experimental results (see Sec. |Vj), and we therefore con- 
clude that including the ^-parameter is not only poorly 
grounded theoretically, but also produces worse quantita- 



tive results than excluding it. For these dots the dot size 
is small enough to have quenched the influence of the re- 
mote bands sufficiently to exclude the remote bands from 
the calculations. 



Appendix B: Averaging methods 

To avoid the complexity of fully quantum mechanical 
calculations, one might consider using simpler averaging 
schemes to calculate the ^-factor. The rationale behind 
such schemes is that an electron or hole experiences a 
spatially varying environment with spatially-dependent 
bulk-like properties. A property £ of this state can be cal- 
culated by weighing the bulk-like property £buik (r) with 
the probability density: 

<£>= [ |*(r)| 2 £ bulk (r)dr 

where \£ (r) is the wave function of the state. Such 
averaging schemes have no validity from a quantum- 
mechanical point of view; perturbing the environment lo- 
cally will alter the whole quantum mechanical state and 
not merely the value of (£). Furthermore, there is not a 
unique choice for £t>uik (r). Some choices might give ac- 
curate results in particular cases, but there is not a single 
choice that would be generally valid. The minimal ap- 
proach to calculate the properties of the electronic state 
properly for these QDs is eight-band k • p-theory. How- 
ever, to stress the differences between the methods, we 
have calculated the electron ^-factor using two different 
averaging schemes. 

For the unstrained electron ^-factors case, we have 
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used an adapted Roth- formula : 



XaX 



#Roth 



2E£A 



3^gap (^gap + A X ) 



(Bl) 



where is the Kane energy and A x is the spin-orbit 
coupling of material X, and E g3iY> the energy gap between 
the ground electron and hole state. The use of E gap takes 
some effects of confinement into account. Penetration 
into the barrier material is taken into account by weight- 
ing with the probability density (r) | 2 : 

<5Roth>= / |*(r)| 2 ^ oth (r)dr, 
Jv 

TnAs TnP 
= WQDfeoth + ^Bfeoth, 

where iuqd,b an d ^Rot^ are the fraction of the proba- 
bility density and the Roth ^-factor inside the InAs QD 
and in the InP barrier material, respectively. Note that 
this information is taken from the k • p wave functions 
at zero magnetic field. The g c -factors obtained via this 
method are indicated by the squares and dotted lines in 
Fig. [71(a). The trend from this averaging method is 
similar to the one from the k • p-calculation. However, 
it underestimates the value of the g c -factor, since orbital 



angular momentum quenching is not taken into account. 
This is even more clearly observed from the discrepancy 
with the quantum wire limits; the lack of orbital angu- 
lar momentum quenching causes the quantum wire limits 
to be much closer to the bulk InAs electron ^-factor of 
- -14.4. 

To show the necessity of including confinement, the 
strained g c -factor case has also been calculated using a 
different averaging scheme. To fully incorporate the ef- 
fects of the strain e (r), one needs to include band mixing. 
The valence bands are mixed due to strain, leading to a 
mixing of the split-off and light-hole states that is linear 
in the strain^: 



|A> 
IB) 
|C> 



3 3 
2' 2 

W- 

2' 2 

±,±± 

2' 2 



ap (e) 
V2 

ap (e) 
V2 



2' 2 
2' 2 



The parameter ao (e) determines the mixture. The g c - 
factor for a strained bulk crystal is then, to first order in 



straii 



>40 



n x / x = o 2E p { 3 1 1 , R2 x 

^ strainedU 3 \2[EcB{e)-E K (e)\ 2 [E C b (e) - E B (e)] E C b (e) - E c (e) 1 } 

1 1 

-a (e) 



E C b (c) - E B (e) E C b (e) - E c (e) 



where Ecb, ^a, ^b, and Eq are respectively the con- 
duction band (CB), the heavy- hole like band (A), the 
light-hole like band (B), and the split-off like band (C) 
edge energies. These energies and the mixing parameter 
ao (e) depend on strain. The g c -factor can then be de- 
termined by weighting the locally varying #^ rained (e (r)) 
by the probability density: 

(^strained) = / |^ 0) ^strained (e (r)) dr. 
JV 



The resulting calculated g c -factor is indicated in 
Fig. [3(b) by the squares and dotted lines. The trend 
is opposite to the k • p-calculation, mainly due to the 
absence of confinement in this averaging scheme. This 
shows that confinement is essential in determining the 
^-factor. 
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